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ABSTRACT 

We report ten lens candidates in the E-CDFS from the gems survey. Nine of the systems are new 
detections and only one of the candidates is a known lens system. For the most promising five systems 
including the known lens system, we present results from preliminary lens mass modelling, which tests 
if the candidates are plausible lens systems. Photometric redshifts of the candidate lens galaxies are 
obtained from the COMBO- 17 galaxy catalog. Stellar masses of the candidate lens galaxies within 
the Einstein radius are obtained by using the z-band luminosity and the V — z color-based stellar 
mass-to-light ratios. As expected, the lensing masses are found to be larger than the stellar masses 
of the candidate lens galaxies. These candidates have similar dark matter fractions as compared to 
lenses in SLACS and COSMOS. They also roughly follow the halo mass-stellar mass relation predicted 
by the subhalo abundance matching technique. One of the candidate lens galaxies qualifies as a LIRG 
and may not be a true lens because the arc-like feature in the system is likely to be an active region 
of star formation in the candidate lens galaxy. Amongst the five best candidates, one is a confirmed 
lens system, one is a likely lens system, two are less likely to be lenses and the status of one of the 
candidates is ambiguous. Spectroscopic follow-up of these systems is still required to confirm lensing 
and/or for more accurate determination of the lens masses and mass density profiles. 
Subject headings: dark matter - gravitational lensing: strong - methods: data analysis - surveys: 
GEMS 



1. INTRODUCTION 

Within the current standard cosmological model of 
the Universe, the gravitationally dominant component of 
matter is dark. Hence, dark matter governs the forma- 
tion and evolution of the luminous matter like galaxies 
in the Universe. Gravitational lensing has proved to be 
a promising technique in the past few decades allowing 
us not only to measure the mass of dark matter halos 
which host luminous galaxies in their centres but also 
to probe the underlying dark matter distribution of the 
lensing halo. 

Lensing has been used, for example, to constrain 
the slopes of density profiles in the central regions o f 
galaxies (e.g., iKoopmans et al.l 120061: iMore et aT1l 200 
and clusters (e.g.. iSand et al.l 120021 : iKneib et al.l 1200 
and to constrain cosmological paramete rs like the Hub- 
ble c o nstant to a better accuracy (e.g.. iFassnacht et al.l 
2002; Coles 2008) and the cosmological constant (e.g., 
Oguri et al.l 12008). Lensing surveys complemented with 
Hubble Space Telescope (HST) observations have ex- 
plored various aspects of galaxy formation and evo- 
lution like studying the interstellar medium of lens 
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galaxies (Fal co et al.l 1 19991) . surface brightness evolu- 
tion (Kochanek et al. 2000), mass density profile evo- 
l ution jK oopman s et al.l 12006). and stellar populations 
(Treu et al. 2006) of early-type galaxies. 

Among the large-area surveys at high spatial reso- 
lution, the GEM S (Galaxy Ev olution from Morpholo- 
gies and SEDs, iRix et aLll2004f ) survey in the Extended 
Chandra Deep Field South (E-CDFS) is the largest con- 
tiguous HST survey with color imaging data to date. 
GEMS was initiated to study the evolution of normal 
and active galaxies out to high redshifts. gems observed 
0.22 sq. deg. inside the E-CDFS in the optical F606W 
(=V) and F850LP (=z) filter band passes with one or- 
bit per pointing. In the center, it contains the GOODS 
south area. Of -60000 objects in the field, -14000 
are galaxies down to R = 24 with photometric redshifts 
from the COMBO- 17 survey. The spectral energy distri- 
butions (SEDs) available from COMBO- 17, besides allow- 
ing accurate photometric redshift estimates, allow clas- 
sifica tion into various galaxy types, including active nu- 
clei ([Wolf et al.l 12004). Stellar mass estimat es for all the 
galaxies are also available from SED fitting (jBorch et al.l 
200|). 

The depth and relatively large volume of the GEMS 
survey is apt for finding incidences of strong gravita- 
tional lens systems. For example, the strong lens survey 
from All- wavelength Extended Groth strip International 
Survey (aegis) which has comparable depth and sur- 
vey area to GEMS, has found a sample of three lenses 
(jMoustakas et al.l 120071 ). Also, the strong lens survey 
from COSMOlogical evolution Survey (cosmos) which 
has a much larger survey area in a sing le band with HST , 
find about 20 good lens candidates (jFaure et~ al. 2009). 
An unbiased sample of lenses can be used to constrain 
the faint end of the background source luminosity func- 
tion, to measure the optical depth for lensing, to study 
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the dark matter fraction in galaxies and its evolution, 
and to study the connection between dark matter halos 
and galaxy properties. 

In this paper, we present a sample of gems lens can- 
didates in the E-CDFS. The structure of the paper is as 
follows. In section 2, we calculate the expected number 
of lenses and explain the steps involved in the selection of 
the lens candidates. In section 3, we describe the decom- 
position of HST images into candidate lensing galaxy and 
candidate lensed images. The decomposition is carried 
out by using the two-dime nsional surface br ightness mod- 
elling program Galfit (|Peng et al.l 12010). We model 
the candidate lens systems as gravitational lenses and 
check whether the derived model parameters are consis- 
tent with a gravitational lens case. If the image configu- 
ration is not due to gravitational lensing, it is likely that 
no sensible lens model can reproduce the lensed image 
configuration. We present results of the mass modelling 
in section 4. In section 5, we compute the stellar mass 
of the candidate lensing galaxy, using color-based stellar 
mass-to-light ratio (MLR) and compare with the lens- 
ing masses as a sanity check. Furthermore, we compare 
the sample properties with those in the literature, and 
discuss the probability of the candidates being lens sys- 
tems. We end with a brief summary in section 6. We 
use the following cosmological parameters in our calcula- 
tions: = 0.3, ^ A = 0.7, and H = 70 km s _1 Mpc -1 . 

2. GEMS LENS SAMPLE 

2.1. Expected Number of Lenses 

The number of lenses expected to be found in a given 
survey depends upon the survey area, the source number 
density, n s , and the optical depth for lensing, r. The 
number of lenses per unit steradian, n\ is given by 



r(z a , ran™) n s (> s , m\ im )dz s , 



(1) 



where z s is the source redshift and m\[ m is the intrinsic 
source limiting magnitude. We estimated the expected 
number of lenses in the GEMS surve y by using optical 
depths from Table 2 of iFaure et all (|2QQ9D . These op- 
tical depths were estimated by integrating the lensing 
cross-section of halos in the Millennium Sim ulation. We 
calcu lated n s by following Eq. 1— Eq. 3 of IFaure et al.l 
(|2QQ9h which requires the distribution of source redshifts 
and source counts. We adopt the source redshift distri- 
bution given by 



p(z 8 ,mn m ) 



$zl 



T(3/f3)z$ 



exp 



(2) 



where f3 = 3/2 and zq = 0.13 mn m — 2.2 and use the cu- 
mulative source counts distribution from the F814W(J) 
band data of the Hu bble Deep Field South (HDF-S, 
iCasertano et al.l 12000). The V band limiting magnitude 
of GEMS is used to calculate a corresponding / band lim- 
iting magnitude. The HDF-S sources follow the relation 
V — 0.49 ~ I a nd GEMS has a V band limiting mag- 
nitude of 28.25 (jRix et al.l [20041 ). Thus, a correspond- 
ing equivalent limiting magnitude in / band for GEMS 
would be 27.76. We note that due to lensing magnifi- 
cation, sources which are intrinsically fainter than the 
limiting magnitude might also fall into our lens sample. 



This would mean that our calculation of expected no. of 
lenses based on sources only brighter than limiting mag- 
nitude might be biased. However, we can de-magnify the 
lensed images, using the magnification factors from the 
lens mass models, to calculate the intrinsic magnitudes 
of respective sources. These source magnitudes happen 
to lie well-above the limiting magnitude used in our cal- 
culation. Therefore, assuming our lens mass models are 
correct, our calculation of expected number of lenses is 
less likely to be biased due to this effect. By integrating 
Eq.Q]from z 8 \ = to z S 2 = 3, we expect about eight lens 
systems in a survey area of 0.22 sq. degrees, down to the 
limiting magnitude. 

2.2. Selection of Lens Candidates 

We adopted the following approach to find gravita- 
tional lens candidates in the GEMS field. As a first 
step, we selected all galaxies above a stellar mass of 
3 x 10 10 and redshift z < 1.1 from the COMBO- 
17 catalog (jBorch et al.l 12006) of 25 000 galaxies. Mas- 
sive lens galaxies are generally more efficient lenses and 
produce larger image separations which are easier to de- 
tect. These galaxies are also, generally, the more lumi- 
nous galaxies at each redshift. As a result, such galax- 
ies have a high signal-to- noise ratio (SNR). The cut at 
z = 1.1 ensures that galaxies with reliable redshifts are 
included in the sample. Also, this limit corresponds to 
the redshift at which the redder z-band filter moves be- 
yond the 4000A-break, so we would start to see pure UV 
images beyond this redshift. In addition to a decrease 
in the SNR, more clumpy star forming structures show 
up more prominently in such cases, increasing the chance 
for contamination of good lens candidates with normal 
star- forming galaxies. These limits helped us to narrow 
down our search to a massive subsample of galaxies which 
are more likely to exhibit strong lensing signatures. The 
first step leaves 1225 COMBO- 17 galaxies of which 1054 
are covered with ACS imaging data by the gems survey. 

In the second step, we visually inspected the 1054 
galaxies, looking at both V — z color composite images 
as well as the two bands individually in different gray- 
scale stretches. This allowed us to see both the surface 
brightness structures as well as color-differences, if any, 
between the image components. We selected possible 
lens candidates based on the presence of a) multiple im- 
age components or arc-like structures mimicking typical 
lensed image configurations, and/or b) color differences 
to the main body of the central galaxy. We found 29 
candidates at the end of this inspection with some be- 
ing more likely lens systems than others. Three people 
ranked these candidates visually from 1 being the least 
likely to 4 being very likely candidates. A total of 19 can- 
didates had average ranks below 2. Of the remaining 10 
candidates, five had average ranks between 2 and 3 and 
five had ranks >= 3. One of the top five candidates is a 
confi rmed lens system (J033238-275653, Blak eslee et al.l 
I2004D with COMBO-17 ID 15422. We analyze this sys- 
tem along with the rest of the candidates. The top five 
candidates with ranks >= 3 are listed in Table [T] and 
examined in the rest of this paper whereas the five low 
probability candidates are described in the appendix. 

3. SURFACE BRIGHTNESS MODELLING 
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We use the publicly available software Galfit to 
model the two-dimensional light distribution in images. 
Galfit is a very flexible code and allows fits to the light 
distribution with a large choice of radial surface bright- 
ness profiles, number of distinct component s, choice of 
PSFs . Galfit in its newest version (V3, iPeng et all 
[20101 ) allows a more flexible analysis through the pos- 
sibility of applying "bending" modes to components and 
adding a banana-shaped bend to an otherwise elliptical 
shape. This facilitates adapting the models of the source 
in the presence of gravitational arcs. 

Generally, we model all systems with one or two com- 
ponents for the lensing galaxy and up to four for the 
lensed images. The choice of component numbers is not 
physically motivated but was adapted to describe the 
whole system well enough so that the candidate lensing 
galaxy and the putative lensed images can be well sepa- 
rated. Testing was done iteratively from repeated mod- 
elling runs with different numbers of components and dif- 
ferent degrees of freedom to each component e.g. bending 
modes switched on or off. The final choice is slightly sub- 
jective but the solution regarding the exact parameters, 
at the same time, is quite degenerate. Since we do not 
use output model parameters from Galfit itself, but 
only from the extracted images, this procedure is very 
robust. We give the individual number of components 
used for each object in Table [TJ 

For four out of five objects, both V and z band images 
are available and we model the deeper V band images 
using the Galfit software. The object 15422 lies on 
GEMS tile #44 for which the V band is missing and thus, 
we revert to modelling the z band. Each row in Fig. [1] 
corresponds to one candidate. From left, the columns 
show a composite V—z color image of the sample object, 
where available, as well as the Galfit modelling results: 
the input V (or z) band image, the surface brightness 
model for the lensing galaxy, and the frame resulting 
from subtraction of the model showing only the lensed 
images. 

4. LENS MASS MODELLING 

4.1. Standard Mass Models 

The magnitude of deflection and distortion of a light 
ray bundle arriving at a point in the lens plane from 
a background source depends upon the surface density 
distribution of mass in the lens. To zeroth order, the 
surface mass density distribution around galaxies follows 
a spherically symmetric singular isothermal profile (SIS) 
given by 



s(0 = 



2£G 



(3) 



where a is the velocity dispersion of the dark matter par- 
ticles, £ is the physical separation in the lens plane and G 
is the gravitational constant. Dark matter halos are tri- 
axial in nature. Therefore an elliptical mass distribution 
is a more realistic description of the surface density dis- 
tribution. The expression that describes the surface den- 
sity distribution for a singular is othermal ellipsoid (SIE ) 
with an axis ratio, q, is given by iKormann et all Jl994), 



■e 2 



(4) 



Here, £i and £2 denote the projected physical separa- 
tions along the major and minor axes of the distribu- 
tion, respectively. The lensing observables are sensi- 
tive to a quantity called convergence (k) which is de- 
fined as the ratio of the surface mass density distribu- 
tion to a critical surface density. The critical surface 
density depends upon the geometrical distances between 
the source, the lens, and the observer and is given by 
S C rit = c 2 D s / (47rG-Dd^ds) 5 where c is the speed of 
light, Dd is the angular diameter distance between the 
observer and the lens, D^s is the distance between the 
lens and the background source, and D s is the distance 
between the observer and the background source. 

Using constraints from extended images is particu- 
larly helpful when the background source is to be re- 
constructed for further analysis. 

We use the software package gravlens (|Keetonll2Q0l[ ) 
to make lens mass models for our candidates. The soft- 
ware package uses the image positions (and optionally, 
their fluxes) and the lens position along with their error 
bars as constraints and outputs a model that best fits 
the constraints. For the purpose of estimating masses 
of the lens galaxies, use of the surface brightness peak 
positions is usually sufficient. A significant improvement 
is not expected in the estimate of the Einstein radius by 
using the extended surface brightness information. This 
is postponed until deeper imaging is obtained especially 
for studying the background sources. 

The convergence at a point (x = ^i/D^, y = ^/Ai) 
in the sky is defined aS 



K(x,y) 



b' 



2gy / x 2 + y 2 /q 2 



(5) 



where b' is the critical radius (also known as the Einstein 
radius). As is clear from Eqs. 2]and[5j the critical radius 
and the velocity dispersion (a) are related by 



(7 




4tt D d 



(6) 



Integrating the surface mass density within the area 
enclosed by the tangential critical curv^EI yields the lens 
mass within the critical curve, Me, which is given by 



M E 



1 



q 2 C 2 D d D s 



2q 4G D d 



(7) 



The first factor on the right hand side of Eq. [7]becomes 
unity in the case of a SIS profile (i.e., when q = 1). In 
addition to the surface mass density distribution, an ex- 
ternal shear is often added to the lens model to account 
for the tidal effects of the environment on the lensed im- 
age observables. The position angles in the lens models 
are measured East of North such that North is deg and 
East is 90 deg. 

An SIS model has a total of 5 free parameters: the 
position (x,y) of the lens, its critical radius, and two 

8 For further details, see eqs. 3.24 - 3.26 in the GRAVLENS 
manual. 

9 The tangential critical curve is an iso-density contour within 
which the (k) = 1. It is the locus of points where lensed images 
are infinitely magnified and are stretched in a tangential direction. 
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N 
▲ 

E ■+ 




15422 

1" 








34244 

• 


o 


• 


m 


40173 

• 








43242 





46446 




Fig. 1. — Five lens candidates of the GEMS survey. The columns are (from left) V-z color composite, V-band image (z-band for object 
15422), surface brightness model of the lens galaxy, and lens galaxy— subtracted residual image showing emission from lensed images only. 
Markers show a scale of \" . The cutouts in the first and the fourth rows are 4.5x4.5 arcsec 2 whereas the rest are 3.0x3.0 arcsec 2 . 
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parameters that describe the position of the source. An 
SIE model has two additional parameters, the ellipticity 
(given by 1 — q) and its position angle. The external shear 
introduces 2 more parameters, the shear strength, and its 
position angle, to each of these models. A doubly imaged 
system provides a total of 6 constraints: 4 constraints 
from the 2 image positions and 2 constraints from the 
lens position. In the absence of any constraints from the 
flux ratios, we only use the SIS model so as to avoid 
an over-fit to the data. Similarly, a quadruply imaged 
system will have a total of 10 constraints: 8 constraints 
from the 4 image positions and remaining 2 constraints 
come from the lens position. This allows us to fit the 
data with SIS+shear and SIE+shear models. In case of 
SIE+shear model, the external shear is often found to be 
degenerate with the ellipticity. 

4.2. Mass Modelling Results 

Using the above models, we present results from pre- 
liminary mass modelling of the lens candidates. The aim 
is to test if simple mass models can reproduce the con- 
figuration of lensed images with reasonable values of the 
model parameters. 

Two kinds of image configurations are commonly seen 
in four-image lens systems, namely, a cusp and a fold. In 
the former, three lensed images are highly magnified and 
tend to merge whereas, in the latter, any of the adjacent 
two images are highly magnified and tend to merge. We 
assume either of the above image configurations (except 
for candidate 40173 which is probably doubly imaged) 
for our mass models and describe the tests carried out 
on all of our candidates. 

Throughout the modelling, we use peak positions of 
the images and the lens galaxy as constraints for the 
mass models. The positions of the images were visually 
determined at the peak of surface brightness (whenever 
possible) whereas the peak position for the lens galaxy 
was chosen by fitting a Sersic profile. The errors on the 
positions of the images and the lens galaxy were assumed 
to be 0.06 arcsec and 0.03 arcsec, respectively. Note that 
the pixel scale is 0.03 arcsec. 

In the following, we describe the mass modelling anal- 
ysis for five of the best candidates. The results of the 
mass models for each candidate are shown with a pair of 
panels in Fig. [21 Every left panel has the model-predicted 
positions overlaid on the true lensed images with the fore- 
ground galaxy subtracted. Every right panel shows a re- 
constructed lensed image configuration assuming a Sersic 
model for the background source with arbitrarily cho- 
sen values for the parameters and using the best-fitting 
model for the lens galaxy. This is for illustrative purposes 
only. The solid green curves in the figure are the critical 
curves corresponding to the best-fitting lens model. 

15422: This is the only system from o ur sample which 
is kn own to be a confirmed lens system (jBlakeslee et al.l 
2004, hereafter, B04). The lens galaxy has a spectro- 
scopic redshift of 0.62 and a photometric redhift of 2.4 
is estimated for the background source. The mass mod- 
elling of B04 suggests that the velocity dispersion of the 
lens galaxy is ~ 305 km s _1 . 

For the sake of consistency and completeness, we make 
mass models with our data and compare our modelling 
results with B04. In this system, one can clearly see 
a thin elongated arc. A faint counter-image is usually 



expected and can actually be seen in Fig.[T]where the lens 
galaxy has been subtracted. This counter-arc is clearly 
visible in B04. We use the peak positions of all four 
images as constraints to the mass models. 

We test an SIE mass model alone motivated by the el- 
liptical light distribution of the lens galaxy. A SIE model 
fits reasonably to this lens system with a reduced x 2 of 
3. The Einstein radius is found to be 1.20 arcsec. The 
mass model suggests an ellipticity of 0.60 with a position 
angle of 64 deg. This is consistent with the ellipticity 
and the position angle measured from the distribution 
of light from the lens galaxy. The image positions pre- 
dicted by the SIE model are shown with circles in the 
left panel of A in Fig. [2] overlaid on the lens-subtracted 
optical image. The right panel of A in the figure shows 
a reconstructed lensed image assuming a source lensed 
by the best-fitting SIE model. The green ellipse shows 
critical curves of this model. 

In conclusion, the SIE model constrained by the arc 
and counter- image provides a reasonably good fit. A 
comparison of our modelling results with B04 is done in 
Sect. El 

34244: In this system, a cusp configuration is likely 
and a counter-image is clearly detected. As before, we 
test three possible mass models namely, SIE, SIS+shear 
and SIE+shear. 

The SIE model suggests an ellipticity of 0.52 with a 
position angle of —31 deg. The Einstein radius is found 
to be 0.40 arcsec. The model predicted image positions 
overlaid on the lens-subtracted residual image are shown 
in left panel of B in Fig. [2] The right panel of B shows 
an artificial source lensed with the best SIE model along 
with the critical curves. A SIS+shear model indicates a 
high shear of about 15 per cent also at the position angle 
of —31 deg and Einstein radius of 0.42 arcsec. 

Both the ellipticity and shear could not be simultane- 
ously constrained with the image positions unless a prior 
on the ellipticity was used from the intensity distribu- 
tion. Thus, to test a SIE+shear model, an ellipticity of 
0.36 with a position angle of —42 deg was used from the 
light profile of the lens galaxy. The best-fit model has a 
moderately high shear of 9 per cent with a position angle 
of —11 deg. 

All of the above models have reduced x 2 < 1.5 and 
hence, are equally favorable. Using the constraints from 
the image positions, the SIE, and SIS+shear models sug- 
gest a position angle of —31 deg for the ellipticity in the 
lens potential. On the other hand, the position angle of 
the intensity distribution of the lens galaxy shows a sig- 
nificant offset from the model predicted value indicating 
inconsistencies in the alignment of matter and light dis- 
tribution. It is unclear if any one of the above models is 
better than the others due to lack of physically motivated 
values of the parameters. 

Here, the SIE model may be referred to as the current 
best-fitting model. A subsequent investigation of the lens 
and its environment might help in testing whether the 
misalignment between the matter and light distribution 
is physical. 

40173: Owing to the circular symmetry of the sur- 
face brightness distribution, this candidate could be an 
example of a partial Einstein ring or a double lens sys- 
tem. Since the lens-galaxy light subtraction may not be 
accurate, the residual arm-like structure could be mis- 
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Fig. 2. — Each candidate is represented by a pair of panels. The left panel in the pair shows the model-predicted positions of lensed 
images and the lens galaxy overlaid on the residual images. The right panel in the pair shows reconstructed lensed images assuming an 
extended source and using the best-fitting lens model. Also, critical curves are shown with green solid curves. 



TABLE 1 

GEMS STRONG LENS CANDIDATES 



ID a 


Tile b 


RA C 


DEC C 


Vtot d 


£tot d 


zf 


^lsp 


comp. g 


log(M t *) h 


\og(M*y 







log(ME) 1 


15422 


44 


03:32:38.21 


-27:56:53.2 




19.20 


0.58 


0.62 


1+3 


11.4 


ii, ^-0.20 


1.20 


0.40 


u on+0.1b 


34244 


94 


03:32:06.45 


-27:47:28.6 


23.43 


21.11 


1.00 


1.02 


1+4 


11.0 


10.82t°;*° 


0.40 


0.48 




40173 


35 


03:33:19.45 


-27:44:50.0 


20.86 


19.67 


0.44 


0.42 


2+3 


10.7 


10.58lg;S 


0.59 


1.00 


io.85l°:S 


43242 


45 


03:31:55.35 


-27:43:23.5 


21.85 


20.00 


0.66 


0.66 


2+3 


11.1 


11.04±g;iS 


0.57 


0.60 


n.08±8:S? 


46446 


47 


03:31:35.94 


-27:41:48.2 


23.60 


21.70 


0.88 




2+3 


10.6 


10.37l°-j3 


0.42* 


1.00 


1 QD+ - 23 



a-ID from the COMBO- 17 catalog. b-GEMS tile number (1-63 GEMS, 80-95 GOODS region), c-positions of candidates are in J200 with 
units of hrs:min:sec and deg:min:sec. d-Apparent V=F606W and z=F814LP HST/ACS magnitudes of the whole system ( Ca ldwell et al. 
2008). e-Photometric redshift of lensing galaxy from COMBO-17. f-Spectroscopic redshift of lensing galaxy (see Sect. 15.21 for details). 
g-Number of components used in G alfit composition for lensing galaxy + lensed images. h-Total stellar mass of the lens galaxy in 
units of Mq from Borch et al. (2006). i- Stellar mass of the lens galaxy inside the critical radius in the units of Mq. j-Einstein radius 
in arcsec corresponding to SIE model except * is corresponding to SIS+shear model, k- axis ratio from the lens mass model 1-Mass 
within the critical curve in units of Mq. The upper and lower limits correspond to the 16 th and 84 th percentile of the lens mass PDF. 



leading. Instead, we choose a two-image configuration. 
An SIS model, centered at the position of the lens galaxy, 
can be tested using the peak positions of two images as 
constraints. An Einstein radius of 0.59 arcsec with a 
symmetric isothermal density profile reproduces the im- 
age separation. Model predicted image positions and an 
ideal source lensed by the mass model are shown in left 
and right panels of C in Fig. [2j respectively. 

However, fitting the image separation alone can not, 
necessarily, indicate gravitational lensing in the system. 
Spectra of the lensed images will be crucial in concluding 
whether blue arc-like feature is a site of star-formation at 
the lens redshift or a lensed arc. Note that, if this indeed 



is a lens system then size of the Einstein radius can give 
a very accurate measure of the mass within the Einstein 
radius, irrespective of the assumed density profile. 

43242: This lens system appears to be in a fold config- 
uration. We assumed that the two southernmost images 
are merging to form an arc- like feature. With the posi- 
tions of four images, following models are tested. 

Firstly, SIE and SIS+shear models are tested. Both 
the models fit equally well with a reduced x 2 ~ 2. The 
Einstein radius of the lens galaxy is about 0.57 arcsec 
with an ellipticity of 0.4 at a position angle of —20 deg 
or alternatively, a shear with 13 per cent strength at the 
same position angle. The panels of D in Fig. [2] show 
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the critical curves for an SIE mass model on the right 
whereas the observed and model-predicted image posi- 
tions (green circles) on the left. The ellipticity of the 
bulge of the lens galaxy is 0.32 with a position angle of 
—25 deg which is consistent with the ellipticity of the 
SIE model. Using this ellipticity as a constraint in a 
SIE+shear model results in equally good fit. The shear 
for this model is 6 per cent with a position angle of 1 deg. 
This implies that the bulge component along with rea- 
sonable shear strength can reproduce the configuration 
of the lensed images. 

Deep observations are certainly desirable to well- 
resolve the arc and high surface brightness images which 
are contaminated due to the emission from the lens 
galaxy. Furthermore, if this system hosts a lens then 
spectroscopic redshifts will provide better constraints in 
measuring the total mass of the late type lensing galaxy 
bulge. This would be interesting since there are only few 
late type lenses known in the literature (e.g.jFeron et all 
120091 iSvgnet et ai1l2Q10t iDutton et aHl2Qllh . 

46446: This candidate indicates a fold configura- 
tion, more obvious in the bottom-right residual image of 
Fig. [U We model the lens system with standard isother- 
mal models. 

The SIS+shear and SIE models for this system suggest 
Einstein radii of 0.42 arcsec and 0.31 arcsec for the lens 
galaxy, respectively. The reduced x 2 f° r these models is 
~ 2. The models need either high shear of 30 per cent at 
a position angle of —76 deg or a high ellipticity of 0.85 at 
a position angle of —79 deg. The brightness profile fitting 
to the lens galaxy also implies an unusually high elliptic- 
ity of ~ 0.95. However, the extremely high shear is un- 
expected for a galaxy-scale lens system. Using priors on 
the ellipticity from the light profile, a SIE+shear model 
is tested. The parameter space of ellipticity and shear 
has a broad degeneracy resulting in a range of equally 
favorable models. 

The left panel of E in Fig. [2] shows the predicted image 
positions with green circles using SIS+shear model over- 
laid on the lens-subtracted image. The right panel of E 
in Fig. [2] shows the critical curves and the reconstructed 
image configuration assuming the best-fitting SIS+shear 
model parameters which well-reproduce the positions of 
the observed lensed images. 

The high ellipticity of the lens from the mass model 
is inconsistent with the morphology and optical proper- 
ties of the central galaxy. The concentration of light in 
the central regions of the lens galaxy does not suggest 
a high concentration of mass which is required for it to 
be a lens. In addition, statistically, we expect more red 
galaxies with star- forming shells than red lens galaxies. 
Therefore, it would be decisive to know the redshifts of 
the blue features. 

5. DISCUSSION 

5.1. Color-based Stellar Masses 

We compute t he st ellar mass M* inside of the critical 
curves (see Sect. 14. 2ft by weighting the luminosity in the 
z-band with the stella r MLR as derived from the (V — z) 
color. According to Gal lazzi fc Belli (j2009[ ). who used 
a large Monte-Carlo library of star formation histories 
to assess the systematic uncertainties in calculating stel- 
lar mass-to-light ratios, a blue rest-frame optical color is 



a good proxy for the MLR. Dust attenuation, metallic- 
ity, and age have degenerate effects on the MLR derived 
from colors and on the luminosity. Incidentally, the ef- 
fects act in opposite directions on MLR and luminosity 
which cancel out to a good extent. As a result, the es- 
timated stellar mass is fairly robust against variations 
in these parameters, of which it is not possible to have 
independent estimates. 

We estimated the stellar masses for our lens galax- 
ies with a Ba y esian statistical approach as outlined in 
iGallazzi et al.l (j2005l ). We created a lib r ary o f model 
SEDs by convolving Bruzual fc Chariot! (j2QQ3h simple 
stellar population (SSP) models of different metallicities, 
ranging from 20 per cent to 2 times the solar metallicity, 
with a Monte Carlo library of complex star formation 
histories (SFH), and dust att enuations. We follo wed the 
prior distribution adopted by iSalim et Itll (120051) and as - 
sumed an initial-mass- function (IMF) of lChabri er (2003). 
The stellar mass of each model is computed by multiply- 
ing the M*/Lf850LP of the redshifted model with the 
observed luminosity Lfssolp of the lens galaxy. Inputs 
for each target are the observed V and z fluxes of the 
candidate lensing galaxy inside the critical curve after 
GALFlT-based subtraction of the lensed images as well as 
the redshift of the lens galaxy. For each galaxy, we built 
the probability density function (PDF) of log (M*/M ) 
by comparing the observed V — z color with the color of 
the models computed at the redshift of the galaxy. Only 
models with formation age younger than the age of the 
Universe, at the redshift of the galaxy, contribute to the 
PDF. The stellar mass of the galaxy is then estimated 
as the median of the PDF with a 68% confidence inter- 
val given by the 16 th — 84 th inter-percentile range. The 
results are listed in Table [TJ 

For the object 15422, we only have the flux at F850LP. 
Hence, no color information can be obtained. Motivated 
by the SO-like morphology of this galaxy, we assume that 
it has a rather smooth SFH and intermediate/old stellar 
population. This translates into selecting models with- 
out a burst in the last 2 Gyr and with 4000A— break 
values typical of SO g alaxies (roughly between 1.4 and 2, 
iPoggianti et a l. 2008), in addition to the constraint that 
the formation age be younger than the age of the Uni- 
verse at the galaxy redshift. To estimate the stellar mass 
of this galaxy, instead of computing the PDF, we sim- 
ply take the median and percentiles of the distribution 
in log (M*/Mq) for the models satisfying these criteria. 

5.2. Lensing Mass Estimates 

The Einstein radii and the axis ratios of the best-fitting 
mass models are listed in Table [T] along with the photo- 
metric and spectroscopic (whenever available) redshifts 
of the lens galaxies. Spectroscopic redshifts of the main 
galaxies we re available for 154 22 (Blakesl ee et al.l [2004) 
and 34244 (jVanzella et al.ll2005l ). Furthermore, as a part 
of Arizona CDFS Environment Survey (ACES) spectro- 
scopic redshifts have been measured for the main galaxies 
in 40173 and 43242 (Cooper et al. 2011, in prep). Thus, 
we use the photometric redshift for 46446 only although 
we note that any flux contamination from the putative 
lensed features might bias the photometric redshift esti- 
mate and hence, the lens mass. Any quantitative esti- 
mation of this bias would require tests with simulations 
which is beyond the scope of this paper. Since the red- 
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shifts of the lensed sources are not known, we use the 
PDF of the redshifts of the background source£3 (see 
Eq. [2] in Sect. 12.11 Thus, masses of the lens galaxies 
within the Einstein radii are calculated by substituting 
above described parameters in Eq. [TJ The median masses 
are given as Me in Table [TJ The error bars reflect the 
16 th and 84 th percentile of the PDF of mass, given the 
lens model parameters and the source redshift PDF. 

Using the Einstein radius, lens and source redshifts of 
B04 in case of 15422, we find that their lens mass within 
the Einstein radius is log(ME) = 11.86 M . This mass is 
consistent with our estimate of log(M E ) = H.80lg;J| M 
from our mass model and assuming a PDF for source 
redshift instead (see Table [TJ. In the case of 40173, we 
note that there is a hint of an inner arc-like feature from 
another source which may or may not be at the same 
redshift as the source lensed into the outer arc. Given 
the quality of the image it is difficult to say which of the 
blue features in the inner and outer arcs correspond to 
which sources. Nevertheless, if we assume that the outer 
arc does have a counter-arc on the inner side, then the 
fractional change in the masses is about 0.1 dex which is 
within the uncertainty arising from the lack of knowledge 
of the source redshift. 

5.3. Comparison of Lensing and Stellar Masses 

Table [T] gives the mass estimates obtained from lens 
mass modelling and stellar masses within the critical 
curve of each system. Comparison of the two masses 
shows that lens mass estimates are consistently higher 
than the respective stellar masses of all the candidates. 
This conforms to the scenario that lensing measures total 
projected mass (dark matter and luminous matter) and 
should be larger than the stellar mass, which is sensitive 
to the luminous mass alone, within the same aperture. 

Early-type galaxies happen to lie on a Fundamental 
Plane (FP) defined by the surface brightness, the effec- 
tive radius, and the velocity dispersion. However, ob- 
servationally the early-type galaxies are found to deviate 
from the naive L oc a 2 R expectation and this deviation 
is referred to as the tilt of the FP. The origin of the tilt is 
thought to be due primarily to the content of dark mat- 
ter in the early- type galaxies (e.g., iPadmanabhan et all 
2004, although variations in stellar MLR and structural 
non-homology can affect the tilt as well, possibly to a 
minor degree). Projected dark matter (DM) fractions 
for lens galax y samples have been calculate d in surveys 
like S LACS (jAuger et al.ll2010D and COSMOS ([Faure et al.l 
|2010[ ). The projected DM fraction is found to be posi- 
tively correlated with the lens mass in early-type lensing 
galaxies as shown in Fig. [3l The projected DM fraction 
(/dm = 1 — M*/Me) is calculated for all the GEMS can- 
didates and shown as a function of the lens mass (Me) 
in Fig. [3j For 15422, we also estimate the /dm using 
the lens mass of B04 and stellar mass from Table [TJ 
This is denoted by label "X" in Fig. [3] and is consistent 
with our estimate within the uncertainties. Comparison 
with the SLACS and COSMOS lenses shows that the GEMS 
candidates roughly follow the same correlation. A more 

10 We note that B04 provide a photometric redshift for the back- 
ground source in 15422. However, we use the source redshift PDF 
as a sanity check of our analysis and compare our lens mass esti- 
mate to that of B04. 
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Fig. 3. — Projected DM fraction as a function of lens mass (Me). 
The projected DM fractions in GEMS candidates are similar to those 
found in the lenses from the SLACS and COSMOS surveys. The label 
"X" shows the location of 15422 using the lens mass from B04. 

conclusive result can be presented when the background 
source redshifts are measured which will accurately de- 
termine the lens masses. 

We further checked where do the candidate lens galax- 
ies lie on the halo mass— stellar mass plane. To test the 
above, we need to compare the mass of the halo with the 
total stellar masses. We calculated the velocity disper- 
sion (cr) using the Einstein radius and Eq. [6] We used the 
o~ to calculate the virial mass (M v i r ) f or each of the lens 
galaxi es following the relation given bv lBrvan fc Normanl 

4 a 3 



M v 



(8) 



Here, H(z) is the redshift dependent Hubble parameter 
and G is the gravitational constant. The virial mass is 
defined to be the mass of the halo within a radius that 
encloses a density contrast of A c with respect to the crit- 
ical density of the Universe. We used the fitting function 
provided bv lBrvan fc Normanl ([ 19981 ) to calculate A c as a 
function of the cosmological parameters and the redshift. 

The halo masses (M v [ r ) of these galaxies thus obtained 
are shown as a function of their stellar masses (M t *) in 
Fig. HI We use total stellar massed from iBorch et al.l 
(2006) which are given in Table [TJ The horizontal er- 
ror bars show an uncertainty of 0.3 dex in the stellar 
mass. The vertical bars represent the masses within the 
16 th and 84 th percentile of the PDF. Note that the error 
bars on the mass do not include any systematic uncer- 
tainties due to the assumption of an isothermal density 
profile or the effects of adiabatic contraction of the halo. 
Observationally, early-type galaxies show evidence for 
the contraction of the halo (see e.g., iSchulz et al.l 2010; 

11 Note that the total stellar masses (Mt*) and stellar masses 
integrated inside the critical curve (M*) are consistent within the 
0.3 dex uncertainties. 
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iDutton e t al. 201( 3) whereas late type galaxies show the 
opposite (jDutton et al.ll2010l ). Note that adiabatic con- 
traction, if present, will cause the density profile to be 
closer to our assumption of an isothermal density profile 
in the inner regions. 

We also plot the halo mass— stellar mass relation in- 
fer red from the subhalo- abundance matching technique 
bv iBehroozi et all (j2Q10h at redshift of (dashed curve) 
and redshift of 1 (solid curve) in Fig. |4] for comparison. 
Their relation at z = is in agreement with results from 
weak lensing ( Mandelba um et "all 120061 ). kinematics of 
satellite galaxies dMore et al.ll2010D, an d studies using 
a SDSS group catalog (jYang et al.ll2009f ) w hile the rela- 
tion a t z = 1 is consistent with results from lDrorv et al.l 
(2009). For 15422, we also estimate the M v i r using the 
lens velocity dispersion of B04 and total stellar mass from 
Table CD This is denoted by label "X" in Fig. H and is 
consistent within the uncertainties. 

It is interesting that within the uncertainties of our 
measurements, all of the five systems roughly follow the 
halo mass— stellar mass relation predicted by the subhalo 
abundance matching technique. Although this does not 
conclusively confirm the lensing nature of these candi- 
dates, it certainly adds more credibility to the candidates 
being lens systems. 

5.4. Cross-identification in Imaging at Other 
Wavelengths 

The HST lens candidates from gems were checked 
for cross-identifications in data at other wavelengths, 
namely X-rays, mid-infrared, and radio. The E-CDFS 
was observed with the Chandra X-ray telescope by 
iLehmer et "all (|2005f ). with Multi-band Imaging Pho- 
tomet er for Spitzer (MIPS) at 24/im bv iPapovich et ah! 
(2004)), and the CDFS was sur veyed in the radio u s ing th e 
Very Large Array (VLA) bv iKellermann et all (|20Q8f ). 
None of our lens candidates were found in X-ray imag- 
ing. However, astrometric comparison with the 24/im 
and radio catalogs yielded one matched candidate. 

Object 34244 from our sample is identified as object 82 
from the VLA sample with an offset of Of! 7 in position. 
The radio object 82 lies at an RA of 03 h 32 m 06 s .45 and 
DEC of -27°47 / 29 // .3 and was measured to have a flux 
density of 80 /iJy at 1.4 GHz. This corresponds to a 
1.4 GHz rest-frame luminosity of 10 23 55 W Hz" 1 . This 
candidate is also detected with MIPS at a flux density of 
315.76 /iJy. 

We inferred the total infrared (TIR) luminosity of 
34244 both from the radio fl ux using the radio-IR corre- 
lation (e.g. JYurTet al. 200 J), and from the 24 fim flux us- 
ing the templates of iCharv fc Elbazl ([200 1) to convert it 
to TIR. Hence, the TIR inferred from either the 1.4 GHz 
radio flux of 80/iJy (Kellermann et al. 2008) or the 24/im 
flux of 315.76 /iJy (jPapovich et al.ll2004fl . is -1-5 xlO 11 
Lq. This is in the range of the TIR luminosity of the 
luminous infrared galaxies (LIRG s). We use the radio- 
SFR and TIR-SFR relations from Eel] (|2003f ) to convert 
the radio and the TIR luminosities into a star-formation 
rate. The star- formation rates derived from the radio 
and IR are 16 M /yr and 80 M /yr, respectively. Thus, 
in accordance with their fluffy and not very sharp ap- 
pearance the blue arc-like features in 34244 could well 
be star-forming shells instead of being lensed images of 
a background galaxy. We note that, recently, an optical 
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Fig. 4. — Halo mass estimated from lensing as a function of stellar 
mass. The dashed and soli d curves show the relat ion at redshift of 
and 1, respectively, from Behroozi et al. (2010). The label "X" 
shows the location of 15422 for the halo mass calculated from the 
lens velocity dispersion of B04 model. 

spectrum of the galaxy was ta ken with VIMOS on th e 
Very Large Telescope (VLT) bv I Silverman et all (j2010f ). 
The slit had a width of 1 arcsec and was oriented along 
north-south. The slit is expected to go through a part of 
the arc-like feature. There is no clear detection of pres- 
ence of a higher redshift object, however, any signal in 
the composite spectrum could have been washed out due 
to the seeing (J. Silverman, private communication). 

5.5. Future Investigations 

In this paper, we have carried out various tests to in- 
vestigate the lensing nature of our candidates. However, 
none of our tests can be considered conclusive and the 
candidates need further investigation. We give directions 
for future work in order to confirm or rule out lensing 
and/or what one could learn from these systems, if they 
are true lenses. 

Although 15422 is a confirmed lens system with a 
known spectroscopic redshift for the lens galaxy, no spec- 
troscopic redshift for the background source is available. 
We note that the lens subtracted images of B04 reveal 
another source at a different redshift which is lensed as 
a quad with a smaller Einstein radius than the arc vis- 
ible in Fig. [U Spectroscopic redshifts of both the back- 
ground sources need to be measured for a detailed and 
sophisticated analysis of this system. Currently, very 
few lens systems are known to have chance occurrence 
of two background source s at different redshi f ts behind 
a ga laxy-scale lens (e.g., iGavazzi et all 120081 : iTu et al.l 
2009). These kind of systems are interesting because 
one can constrain the slope of the mass density pro- 
file and cosmologica l parameters to a better accuracy 
(jGavazzi et al.ll2QQ8D . 

The SFR calculations of the main galaxy in 34244 sug- 
gest that it is a LIRG and the arc-like features could be 
sites of star formation. The orientation of ellipticity from 
our mass models are misaligned to the light distribution. 
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TABLE 2 

Comparison of characteristics of recent lens surveys 



Name 


Reference 


Telescope 


Bands 


Area 


Depth 


lens 


Confirmed 










(deg 2 ) 


mag 


candidates 


/followed-up systems 


AEGIS 


Moustakas et al. (2006, 2007) 


HST 


V, I 


0.18 


V = 28.75 


7 


3/0 


GEMS 


This work 


HST 


V, z 


0.22 


V = 28.25 


10 


1/0 


COSMOS 


Faure et al. (2010) 


HST 


I 


1.80 


I = 25.00 


88 


4/18 


CFHTLS-sl2s 


Gavazzi et al., in prep 


CFHT 


u,g,r,i,z 


170 


g = 25.47 


330 


40/65 


SDSS-SLACS 


Bolton et al. (2006) 


Sloan 


spectroscopy 


3732 


r 17.77 


131 


63/131 



Although such misalignments have been noted earlier, 
this might also imply that 34244 is less likely to be a 
lens system. Furthermore, an optical spectrum of the 
candidate suggests no significant detection of a higher 
redshift object which needs to be confirmed with a high 
resolution — high SNR spectrum. 

The current imaging and mass modelling are insuffi- 
cient to comment on the likelihood of lensing in 40173. 
It is not clear if the arc-like feature corresponds to a back- 
ground galaxy unless it is confirmed by spectroscopy. 

The candidate 43242 shows two lines-of-sight to the 
same background galaxy wherein one of the images can 
be seen rather unobstructed and the others are seen 
through the dust in the foreground galaxy. Such lensed 
images can be used to estimate the amount of differ- 
ential dust extinction from the foreground lens galaxy 
(Ehasdot tir et al.ll2QQ6[ ). This system could also be used 
to constrain the disk mass in the lens galaxy which is a 
near edge-on spiral. 

Although the model-predicted positions of lensed im- 
ages in 46446 can be reproduced and they show similar 
colors, the morphology is complex and perhaps contami- 
nated by dust in the lens galaxy. Spectroscopy and obser- 
vations at other wavelengths of this system can shed some 
light on the spectral type of the candidate lens galaxy 
and confirm if it is a star- forming galaxy or not. 

5.6. Comparison with other surveys 

In this section, we briefly compare the efficiency of de- 
tection of lens systems in gems with that of the detection 
of lenses in other surveys. Table [2] summarizes the char- 
acteristics of the surveys like AEGIS, COSMOS, CFHTLS- 
SL2s, SDSS-SLACS along with the gems survey. Lens 
candidates in AEGIS, COSMOS and gems were selected 
from photometric data from HST while the candidates in 
CFHTLS-SL2S were selected from photometric data from 
ground-based telescope. SDSS-SLACS on the other hand, 
used spectral information of galaxies from the SDSS and 
selected candidates based on the presence of emission 
lines at redshifts higher than that of the main galaxy. 
Another sample of ten lens can didates was reported from 
the HST Medium Deep Survey ([Ratnatunga et al.lfl999l ). 
However, we do not include their results in our compar- 
ison because their sample is compiled from a set of 400 
random fields observed with varying limiting magnitudes 
which complicates the selection function. 

The number of lenses found in ground-based surveys 
is large thanks to their large areal coverage. However, 
the small area, high resolution, deep surveys carried out 
with the HST are able to find a large number of lens 
candidates per square degree. The exquisite resolution of 
these surveys allows them to find lens systems with small 
image separation. These lens systems probe the lower 
end of the halo mass function as opposed to the average 



lenses found in low resolution wide-imaging surveys from 
the ground. Such systems are often missed by ground- 
based surveys due to de-blending issues unless one uses 
spectral information, as was done by the SLACS team. 
The depth of the HST surveys also allows detection of 
fainter lensed images around galaxies and also, lenses 
that lie on average at higher redshifts. 

Among the surveys carried out with the HST, COSMOS 
has the largest area, however it has imaging only in a sin- 
gle band, gems and AEGIS have a smaller areal coverage 
but have imaging information in two distinct bands. The 
multi-band information allows a more robust detection of 
lensed images with distinct colors with respect to the lens 
galaxy. In this respect, the gems and the AEGIS samples 
are expected to have a better completeness compared to 
COSMOS. Also, the imaging for COSMOS was carried out 
in the / band which is not particularly suited for detect- 
ing lensed images given that they are usually fainter in 
redder bands. On the other hand, in terms of the fidelity 
of the samples, multi-band information is a double-edged 
sword. As can be seen from the system 34244 in this pa- 
per, shells of star formation in the candidate lens galaxy 
could also be confused to be lensed images thus affecting 
the purity of the sample. 

The AEGIS survey has a very similar survey design in 
terms of area, depth and wavelengths of imaging com- 
pared to gems. Unless the surveys are severely affected 
by cosmic variance, we expect to roughly find similar 
number of lenses in these surveys. The visually identi- 
fied candidate list from AEGIS was further matched to the 
spectroscopic information from the DEEP2 survey which 
reduced false positives in their case. 

We note that a quantitative comparison between the 
completeness, efficiency and fidelity of lens samples in 
different surveys is difficult, primarily due to inherent dif- 
ferences in the design of the surveys and methods used to 
select the candidates. Large uncertainties in our current 
understanding of the faint end of the galaxy luminosity 
function and the redshift distribution of sources compli- 
cate the quantitative comparison further. 

6. SUMMARY 

The GEMS survey consists of HST imaging in V 
and z bands within the E-CDFS covering an area of 
0.22 sq. deg. We searched for lens candi dates from a sam - 
ple of 25 000 galaxies in COMBO- 17 (|Borch et al.ll2006f ). 
After applying our initial selection criteria, a total of 
1054 candidates were selected for visual inspection and 
ranking. This resulted in ten lens candidates (that is, 
nine new candidates and one known lens system) from 
the gems survey, five of which are best cases and were 
further investigated. 

The lensing masses of these five candidates were cal- 
culated under the assumptions of standard density pro- 
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TABLE 3 
Lower rank gems candidates 



ID 


RA 


Dec 


Tile 


V 


z 


^phot 


12589 


03:31:24.89 


-27:58:07.0 


17 


23.36 


21.68 


0.79 


43797 


03:31:31.74 


-27:43:00.8 


47 


23.72 


21.94 


1.02 


28294 


03:31:50.54 


-27:50:28.4 


33 


22.50 


20.32 


0.66 


36857 


03:31:53.24 


-27:46:18.9 


38 


22.19 


20.44 


0.42 


36714 


03:32:59.78 


-27:46:26.4 


37 


20.92 


19.73 


0.44 



files and a PDF for source redshifts. We estimated stel- 
lar masses within the Einstein radii found from the lens 
models. The stellar masses are consistently lower than 
or equal to the lensing masses, as expected. We com- 
pared the projected dark matter fractions from galaxy- 
scale lenses from surveys like SLACS and COSMOS to those 
in the gems lens candidates. Although uncertainty in the 
background source redshifts leads to weak constraints 
on the estimated parameters, we find that all of the 
candidates roughly follow the trend shown by lenses in 
SLACS and COSMOS. Also, these candidates follow the 
halo mass— stellar mass relation of Be hroozi et al.l ([2010) 
within the uncertainties. 

15422 is the only known and confirmed lens system in 
our sample. Our mass estimate is consistent with that 
from B04. However, the imaging in B04 indicates pres- 
ence of a second lensed source in 15422 which deserves 
a detailed analysis. We note that the main galaxy of 
34244 is detected in radio, and the TIR luminosity sug- 
gests that it is a LIRG and that the blue arc-like fea- 
tures might be star- forming shells. The model ellipticity 
of the main galaxy of 46446 is unusually high and in- 
consistent with the light distribution. Results from the 
analyses of the candidates presented here suggest that 
43242 is very likely lens system, 46446 and 34244 are less 



likely to be lens systems whereas the nature of 40173 is 
not clear. Spectroscopic observations of the candidates 
are ultimately needed to definitively confirm or exclude 
lensing in these systems. Also, redshifts of the back- 
ground galaxies will allow the determination of the mass 
of the lens more correctly. Our calculations imply that 
there should be about eight lens systems in gems. Given 
the large uncertainties in the assumptions made in these 
calculations, it is broadly consistent with the two good 
systems or ten plausible candidates in the sample. 
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APPENDIX 
FIVE LOW PROBABILITY CANDIDATES 

The five candidates with average ranks from 2 to 3 are presented here. Color images of these candidates are shown in 
Fig. EJ The candidates are indicated with their identification number from COMBO- 17. Further information like their 
positions, tile numbers, and photometric redshifts of the candidates is given in Table 03 As before, the identification 
numbers and the photometric redshifts are from COMBO- 17 cat alog. All candidates a re on the gems tiles which range 
from 1 to 63. The V and z— band magnitudes in AB are from Cald well et al.l (j2008l ). 
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roughly 6x6 arcsec 2 . North is up and East is left. 
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